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ADDENDUM  -  January  1,  1959 

This  report  was  originally  issued  as  a  classified  docxiraent 
on  December  1,  1957«   No  changes  have  been  made  in  the  manu- 
script, but  some  minor  corrections  and  additional  references 
are  given  below. 

read  "35"  for  "l5" 

"Q  =  10^^  RH" 

read  "Pig.  Ij-c"  for  Pig.  5c" 

read  "Pig,  I+d"  for  "Pig.  5d" 

"  r  =  10'^  hrVt 

read  "Por  example,  with  H  =  10,000  gauss, 
R  =  10  cm,  T  =  10  ev,  we  have   T* =  0.01  sec. 
Or  with  H  =  100,000  gauss,  R  =  100  cm, 
T  =  10  kev,  we  have  V  =   0.01  sec.  again" 

p.  33    eq.(6):     "P^  =  0.08  RTH" 

p.  3I4-  Table  is  somewhat  modified  according  to  the 

above  changes. 

Section  7:         Since  this  report  was  issued,  high  density 

plasmas  (for  which  cusped  containment  is 
appropriate)  have  been  achieved  experimentally; 
e.g.   A.  C.  Kolb,  Magnetic  Compression  of 
Shock  Preheated  Deuterium,  Paper  submitted 
to  the  second  United  Nations  International 
Conference  on  the  Peaceful  Uses  of  Atomic 
Energy. 

p.  1|2  1.  5  from  bottom  and  following:   The  monocusp  is  the 

optimum  shape. 
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PREFACE 

This  report  Is  a  sxirvey  of  and  introduction  to  some  of  the 
more  important  theoretical  aspects  of  those  cusped  geometries 
which  may  be  feasible  as  thennonuclear  devices.  The  possibilities 
are  still  far  from  being  completely  explored,  and  the  situation 
may  still  be  expected  to  change.   Quantitative  analyses  of  many 
of  the  specific  problems  mentioned  in  this  report  will  follow, 
but  in  view  of  the  large  number  of  inter-related  phenomena  which 

must  be  simultaneously  analyzed,  these  results  should  be  considered 
to  be  tentative  xintll  verified  experimentally. 

The  text  of  this  report  is  essentially  identical  to  a  draft 
which  was  written  in  October  1956.   It  was  felt  to  be  advisable 
to  issue  this  material  in  a  tentative  form  rather  than  to  delay, 
possibly  indefinitely,  for  the  completion  and  improvement  of 
the  theory. 
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1 .   Introduction 

The  point  of  departure  for  these  Investigations  Is  a  certain 

Idealized  plasma -field  configuration.  This  Is  a  free  boundary 

surface  or  Interface  separating  a  plasma  domain  In  which  there 

Is  no  electromagnetic  field  from 

a  vacuvun  domain  In  which  there  Is 

no  plasma.  The  original  advantage 
Vacuum 

of  this  model  was  Its  mathematical 


simplicity  which  allows  deeper  study 
Plasma    \  ° 

and  more  quantitative  analysis  than 
Is  easily  possible  when  plasma  and 
field  Intermingle.   Subsequently, 
Pig.  1  however.  It  was  found  that  such  a 

sharp  separation  (provided  that  Is  can  be  realized  experimentally) 
may  have  definite  practical  advantages  over  similar  configurations 
without  the  sharp  boundary  {  Ref.  [12]). 

Among  the  theoretical  questions  that  should  be  asked  about 
a  given  configuration  presumed  to  be  operating  as  a  thermonuclear 
device  are  Its  stability,  rate  of  loss  of  particles  and  of  energy, 
and  experimental  feasibility.   In  order  to  answer  all  but  the 
first  question,  the  original  Idealization  must  be  modified  to 
take  Into  account  a  continuous  transition  from  plasma  to  vacuiim. 
However,  the  plctiore  we  have  In  mind  la  still  that  of  a  thin 
transition  zone  of  mlnlmxira  theoretical  thickness  separating  a 
reasonable  approximation  to  a  vacuum  from  an  almost  field-free 
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plasma.   It  Is  this  feature  that  distinguishes  the  configurations 
analyzed  in  this  report.  Treatment  of  the  general  case  will  be 
left  to  later  reports  but  will  be  occasionally  referred  to  for 
comparison. 
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2.  Description  of  the  Geometries 

The  basic  two-dimensional  stable  geometries  are  Indicated 
In  Pig.  2.  These  shapes  have  been  computed  analytically  by  con- 
formal  mapping  (Ref.[l4]).  The  essential  mathematical  feactire  of 


(b) 


Fig.  2 
the  pz*oblem  of  the  determination  of  such  Interfaces  Is  that  the 
position  of  the  Interface  Itself  must  be  found  as  well  as  the 
accompanying  magnetic  field  configuration,  but  the  latter  Is 
determined  by  the  surface  currents  on  the  unknown  Interface  as 
well  as  by  the  given  external  colls.   An  equlllbrlxim  free- 
surface  configuration  Is  distinguished  from  any  other  magnetic 
field  surrounding  a  perfect  conductor  by  the  fact  that  the  magni- 
tude of  B  must  be  a  constant  on  the  entire  surface  of  the 
plasma.  This  restricts  the  possible  shapes.  The  special  case 
of  Fig.  2a  In  which  the  external  colls  (Infinite  line  currents) 
are  placed  at  the  comers  of  a  square  and  move  out  to  Infinity 
yields  a  hypocyclold  as  the  free  boundary,  :^'^   +  y^^-^  »  1  . 
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In  actual  practice,  the  colls  will  have  finite  dimensions,  and 
numerical  computation  will  have  to  be  resorted  to.  This  type  of 
computation  is  feasible. 

For  a  given  external  coil  system,  there  is  a  one  parameter 
family  of  equilibrium  free-surface  configurations.   One  can 
imagine  that  plasma  is 'poured  In',  increasing  the  size  of  the 


(a) 


(b) 


Pig.  3 
diamond  xintll  a  certain  maximum  size  is  reached  (at  which  point 
the  cusps  extend  to  slightly  beyond  the  coil  position.  Fig.  Jib), 
after  which  the  plasma  'spills  out',  and  there  is  no  bounded 
solution,  in  other  words,  no  containment. 

The  shape  of  the  interface  is  unchanged  If  we  constiruct  a 
crossed-field  configuration  by  inserting  in  the  plasma  a  constant 
longitudinal  magnetic  field  which  is  perpendicular  to  the  plane 
of  the  external  field  (Fig.  4a);  this  field  is  supposed  to  be 
confined  to  the  plasma  and  does  not  extend  to  the  vacuum  domain. 

In  three  dimensions,  there  are  a  n\imber  of  cusped  configurations 
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which  are  suggested  by  Fig.  2.   We  say  'suggested'  because  the 
conformal  mapping  technique  Is  not  available  In  three  dimensions 
and  the  existence  of  these  configurations,  although  plausible,  is 
not  yet  proved.   In  some  cases,  the  existence  has  been  empirically 
demonstrated  by  the  apparent  convergence  of  numerical  computation 
(Ref.  [15]). 

A  number  of  the  three-dimensional  possibilities  are  sketched 
In  Fig.  4.   At  (a)  we  have  merely  reproduced  Pig.  2a  In  three- 
dimensional  perspective.   The  external  coil  system  consists  of 
four  parallel  wires.  This  can  be  compared  with  a  conventional 
linear  pinch  or  with  a  B  -pinch  if  the  indicated  crossed  field 
is  inserted  in  the  plasma.  We  obtain  (b)  by  bending  (a)  to  close 
the  figure  into  a  torus.   The  shape  of  the  cross-section  will  be 
somewhat  modified  by  the  bending  and  more  so  if  there  is  an 
interior  crossed  field.  The  external  coil  system  is  now  four 
circular  wires.   It  is  interesting  to  note  that  if  we  invert 
the  degrees  of  importance  which  we  attach  to  the  external  cusp 
field  and  to  the  internal  crossed  field,  i.e.  consider  the 
internal  field  to  be  basic  and  the  cusped  field  auxiliary,  then 
this  configuration  becomes  similar  to  some  of  the  helical 
Stellarator  geometries.   Figure  4c  is  2a  rotated  about  an  axis 
through  the  cusps;  this  is  best  compared  to  the  Mirror  Machine. 
The  external  current  system  could  be  two  opposed  Helmholtz  coils. 
Figure  4d  is  Figure  2b  rotated  about  a  horizontal  axis  or, 
equivalently,  a  number  of  4c 's  placed  end  to  end.   The  coll 
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Crossed 
field 


(a) 


(c) 


(e) 


(b) 


(d) 


(f) 


(g) 


Pig.    4 
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(h) 


system  Is  a  succession  of  closed  loops  encircling  the  plasma. 
In  4e  we  bend  4d  around  to  get  a  closed  configuration.  These 

have  been  called  hexacusp,  octocusp,  etc.  depending  on  the  num- 

♦ 
ber  of  segments,  and,  by  analogy,  4c  Is  a  monocusp. 

Figures  4a  and  4d  can  be  fxirther  modified  by  cutting  out  a 
hole  In  the  plasma  using  an  Internal  conductor  as  Indicated  In 
4f  and  4g  respectively.  The  purpose  of  this  Is  economy  of 
magnetic  energy  and  will  be  discussed  later.  This  modification 
also  provides  a  means  of  Introducing  a  crossed  field  Into  4d 
or  4g  by  allowing  the  field  of  this  central  conductor  to  pene- 
trate somewhat  Into  the  plasma.  In  Fig.  4h,  we  have  sketched 
such  a  modification  of  kc   with  a  crossed  field.  The  possibility 
of  cutting  a  hole  In  a  torus  such  as  4b  or  4e  depends  on  passing 
conductor  leads  Inside;  this  also  seems  to  be  theoretically 
possible,  but  the  sketch  Is  rather  formidable.  The  more  complex 
of  the  above  suggestions  are,  of  course,  only  tentative,  subject 
to  either  theoretical  or  experimental  confirmation. 

N.  Levlne  has  proposed  certain  configurations  which  bear 
a  resemblance  to  that  of  Fig.  2a,  but  with  six  colls  Instead  of 
four. 

All  of  these  geometries  are  Idealized  In  that  a  sharp 

separation  surface  Is  postulated.  There  Is  a  corresponding 

family  of  more  general  configurations  exhibiting  a  continuous 

* 
This  was  originally  referred  to  as  the  cuspidor,  reference  [12]. 
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transition.   For  example,  with  plasma  and  field  Intermingled, 
Fig.  2b  represents  the  original  Picket  Fence  configuration 
(Ref.  [11]). 
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3.   stability 

There  Is  one  Important  feature  shared  by  vii*tually  all  the 

geometries  Just  mentioned,  and  that  is  absolute  stability  against 

« 
all  perturbations.   In  the  absence  of  a  crossed  field,  this 

property  has  been  proved  using  the  energy  principle  described  in 
Ref.  [1]  and  with  crossed  fields  using  the  more  general  energy 
principle  of  Frieman  (Ref.  [J>]).     The  crux  of  the  argxjment  which 
demonstrates  this  stability  is,  roughly  speaking,  that  the  mag- 
netic potential  (or  some  equivalent  potential)  is  a  minimum  for 
the  configiirations  of  Section  2.  With  only  slight  modifications, 
the  argument  can  be  applied  to  any  conventional  conducting 
fluid  and  to  any  conducting  gas,  e.g.  one  governed  by  the  con- 
ventional equations  of  gas  dynamics  or  by  the  Boltzmann  equation. 
By  suitable  modifications  of  this  argument,  stability  can  be 
demonstrated  to  hold  for  finite  pertxirbations  as  well  as  for 
infinitesimal  and  for  a  variety  of  external  constraints  including 
certain  external  circuitry  (Ref.  [7]).   The  sole  property  that 
the  fluid  (or  collection  of  particles)  must  possess  is  that  the 
field  and  fluid  do  not  mix;  i.e.  the  transition  zone  must  remain 
thin  for  an  appreciable  length  of  time.   In  particular,  it  should 
be  emphasized  that  the  cusp  itself  is  stable  against  displacements 
or  deformations. 

All  the  geometries  of  Section  2  are  stable  whether  or  not 
there  are  additional  external  fixed  conducting  shells  or  internal 


'References  [2],  [4],  [7]. 
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crossed  fields  but  not  necessarily  when  accelerated.   Inward 
(compressive)  acceleration  tends  to  destabilize  and  outward 
(expansive)  stabilizes.   In  a  periodically  oscillating  plasma, 
the  acceleration  will  be  stabilizing  when  the  plasma  is  most  com- 
pressed and  destabilizing  when  it  is  expanded.   In  the  absence 
of  a  crossed  field,  the  maximum  inward  acceleration  that  can  be 
sustained  is  given  by 


(1)  2aR  <  k^   =  B-^/np 


where  a  is  the  acceleration,  R  is  the  maximiam  radius  of 
curvatxare  of  a  magnetic  line  on  the  interface,  and  A   is  the 
Alfven  speed  constructed  from  the  plasma  density  p  and  vacuum 
field  value  B^  at  the  interface  (see  [^]).  The  maximum  value 
of  R  is  fovind  midway  between  cusps  and  is,  generally,  about 

equal  to  the  radius  of  the  plasma. 

« 
Crossed  fields  have  an  additional  stabilizing  effect  ,  and 

this  effect  is  enhanced  by  surrounding  the  plasma  with  a  conducting 
shell;  but  without  the  crossed  field,  the  external  conductor  has 
no  effect  at  all  on  stability.  The  relative  magnitudes  of  these 
effects  on  stability  can  be  described  by  comparing  them  as  equiva- 
lent accelerations.   In  the  presence  of  crossed  fields,  the 
results  vary  with  the  geometry,  and  we  shall  not  attempt  to  cover 


This  is  rigorously  true  if  the  crossed  field  is  uniform,  e.g.  in 
Fig.  4a,  but  is  an  oversimplification  if  the  crossed  field  is 
curved  as  in  Pig.  4b. 
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all  cases.   For  example,  In  the  geometry  of  Pig.  4a,  the  crossed 
field  effect  has  order  of  magnitude 


(2)  a^-  A^Vl 


where  A^  is  the  Alfven  speed  associated  with  the  internal 
field,  B-j^,  and  L  has  the  order  of  the  length  of  the  plasma 
This  should  be  compared  with  the  curvature  effect  which  is 


(3)  a^  -  A//2R 


The  maxlmvim  acceleration  that  a  given  config\iratlon  can 
withstand  and  remain  stable  is  given  by  the  sum  a.  +  a   .   Of 
course,  in  an  actual  experiment,  the  duration  of  a  destabilizing 
acceleration  would  also  be  a  significant  parameter. 
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4.   Particle  Losses 

What  we  wish  to  discuss  In  this  section  Is  the  containment 
time  of  a  typical  stable  cusped  geometry.  The  perfectly  con- 
ducting fluid  model  with  a  mathematically  sharp  Interface  has 
an  Infinite  time  of  containment  since  It  Is  stable  and  there  Is 
no  diffusion  across  the  field.  This  fluid  model  could  be  modi- 
fled  macroscoplcally  by  Introducing  finite  conductivity,  but  it 
is  more  appropriate  to  do  a  direct  particle  analysis. 

First,  consider  those  cases  in  which  there  is  no  crossed 
field.  The  diffusion  of  particles  across  the  field  turns  out  to 
be  a  negligible  loss  compared  to  the  lcs3  of  particles  through 
the  cusps  in  the  direction  of  B.  The  latter  loss  is  a  direct 
consequence  of  adopting  a  finite  transition  layer  rather  than 
an  abrupt  discontinuity  between  plasma  and  field.   It  is  not 
immediately  related  to  a  mirror-type  loss  since  the  magnitude 
of  B  does  not  increase  for  a  particle  approaching  a  cusp; 
however,  in  a  certain  generalized  sense,  the  mechanism  is  the 
same  as  in  a  mirror  (ftef.  [13]).   The  magnitude  of  this  loss  is 
fovmd  to  be  directly  related  to  the  thickness  of  the  transition 
layer,  and  the  best  result  is  obtained  with  the  thinnest  possible 
layer.   The  gyro  radius  represents,  roughly  speaking,  the  shortest 
distance  within  which  a  particle  can  be  turned  back  by  a  magnetic 
field;  this  represents,  therefore,  a  lower  bound  on  the  thickness 
of  the  layer.   There  is  still  a  considerable  range,  however, 
between  the  electron  and  the  ion  gyro  radii  (a  factor  of  60  for 
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deuterium).   It  can  be  made  plausible  by  appropriate  theoretical 
arguments  that,  by  Judicious  experimental  control,  the  minimum 
theoretical  thickness,  i.e.  the  electron  gyro  radius,  can  be 

approximated  quite  closely.   It  is  this  thin  transition  zone 

* 
which  seems  to  give  the  high  density  (large  p  )  cusped  con- 
figurations an  advantage  over  the  low  density  picket  fence 
conf  igxu*ations . 

The  estimation  of  particle  losses  requires  (among  other 
things)  the  solution  of  an  extremely  complex  self -consistent 
field  problem  involving  the  determination  of  the  electromagnetic 
field  in  the  boundary  layer  (including  the  neighborhood  of  the 
cusps)  30  as  to  be  consistent  with  the  particle  trajectories 
which  give  rise  to  this  electromagnetic  field.  What  we  shall 
do  is  present  a  sequence  of  models  of  increasing  complexity 
until  we  obtain  one  which  we  believe  gives  a  reasonable  approxi- 
mation to  the  actual  loss  rate.   The  details  will  be  presented 
in  later  reports. 

For  the  time  being,  to  simplify  the  arguments,  we  assume 
that  ions  and  electrons  have  the  same  temperature  and  a  mono- 
energetic  (i.e.  fixed  speed)  isotropic  velocity  distribution. 
Any  plasma  which  is  dense  enough  to  be  interesting  thermonuclear- 
wise  must  be  neutral  on  the  whole;  this  does  not  eliminate  the 
possibility  of  non-neutrality  in  very  thin  regions  on  the  order 


* 
P  is  the  ratio  of  plasma  energy  to  magnetic  energy. 
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of  the  Debye  length,  but  docs  require  neutrality  in  any  domain 
as  large  as  the  gyro  radius,  e.g.  in  the  bo\indary  layer. 

The  first  step  Is  to  take  the  magnetic  field  configura- 
tion of  the  fluid  model  with  a  dlscontlnxilty  In  B  and  compute 
particle  trajectories,  cf.  Fig.  5a.  The  trajectories  are 
straight  lines  Joined  by  cycloldal  pieces.  The  limiting  case 


(a) 


(b) 


Pig.  5 

of  zero  gyro  radius  Is  a  billiard  ball  model,  cf.  Pig.  5b.   In 

this  Idealized  case,  every  particle  will  be  turned  back  before 

reaching  the  cusp  unless  it  is  on  the  axis  and  aimed  directly 

at  the  cusp.  There  is  no  loss  of  particles;  the  containment  time 

is  infinite  which  is  consistent  with  the  original  fluid  model.   On  the 
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other  hand,  the  'soft*  billiard  table.  Fig.  3a,  does  allow 
certain  particles  to  reach  the  cusp  and  escape;  (this  will  be 
adopted,  temporarily,  as  the  definition  of  'escape').  More 
precisely,  those  particles  which  are  aimed  near  enovtgh  to  the 
cusp  are  lost.   Analysis  of  these  trajectories  leads  to  the 
conclusion  that  the  rate  of  loss  of  particles  of  a  given  speed 
is  proportional  to  the  mean  distance  of  travel  toward  the  cusp 
in  a  single  encounter  with  the  boundary  layer  ("d' in  Pig.  5a)} 
this  however,  is  proportional  to  the  gyro  radius.  Since  elec- 
trons are  faster  by  a  factor  60  and  have  a  smaller  gyro  radius 
by  the  same  factor,  the  loss  rates  of  ions  and  of  electrons  az^ 
equal  on  the  basis  of  this  conqsutation. 

This  computation,  using  the  original  discontinuous  magnetic 
field,  is  not  self -consistent  in  two  ways;  the  finite  penetration 
distance  produces  a  current  layer  rather  than  a  current  sheet, 
ultimately  yielding  a  continuous  transition  of  the  magnetic 
field,  and  the  greater  penetration  of  ions  produces  an  electric 
field  which  further  modifies  the  trajectories.  Of  these  two 
effects,  the  first  produces  only  minor  quantitative  changes,  but 
the  second  is  crucial.  The  reason  for  this  is  that  charge 
separation  of  an  amount  indicated  by  the  difference  between 
electron  and  ion  gyro  radii  woxild  produce  gigantic  electric 
fields.  To  a  vez*y  good  approximation,  the  electric  field  must 
be  such  as  to  produce  equal  mean  penetration  of  ions  and  of 
electrons,  or  more  precisely,  equal  densities  of  both  particles 
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evex*ywhere  in  the  layer.   Solution  of  this  self-consistent  problem 
yields  the  result  that  the  electrons  are  pulled  out  only  slightly 
but  the  ions  are  pulled  in  to  almost  the  original  electron 
penetration  distance.  Although  the  mean  penetrations  are  now  the 
same  for  the  two  particles,  the  mean  distances  travelled  in  the 
direction  of  the  cusp  are  not  quite  the  same  since  the  actual 
trajectories  of  ions  and  electrons  are  no  longer  similar  in  the 
presence  of  the  electric  field.  However,  a  very  rough  estimate  of 
the  loss  rates  can  be  computed  by  assuming  that  electron  and  ion 
orbits  are  geometrically  similar.  Since  the  electron  penetration 
is  roughly  the  same  as  in  the  original  non-self-consistent  model, 
while  the  ions  have  been  pulled  in  by  a  factor  6o,  we  conclude 
that  the  electron  loss  rate  is  roxighly  the  same  as  was  originally 
computed,  but  the  ion  loss  rate  has  been  cut  down  60  times  by  the 
charge  separation  electric  field.  More  precise  computation 
yields  the  factor  1^  rather  than  60.   (These  computations  will 
be  made  in  succeeding  reports). 

This  second  model  is  internally  self-consistent,  but  the 
difference  in  electron  and  ion  loss  rates  must  be  compensated 
externally  by  a  change  in  the  potential  of  the  whole  plasma. 
An  equilibrium  will  be  reached  in  which  the  ion  and  electz*on 
loss  rates  are  eq\ial,  but  the  crucial  question  is  whether  this 
is  accomplished  mostly  by  increasing  the  ion  loss  rate  or  by 
decreasing  the  electron  loss  rate.  The  answer  to  this  question 
involves  solving  another  self-consistent  problem  involving  the 
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above-mentioned  boxmdary  layer,  the  cusp  neighborhood  where 
two  boundary  layers  merge  and  an  external  (possibly  penetrating 
the  plasma)  field  necessitated  by  the  xonequal  loss  rates  com- 
puted in  its  absence.  The  basic  issue  is  as  follows.  The  ex- 
ternal electrostatic  potential  accelerates  the  emerging  ions 
to  the  walls  but  it  turns  back  a  certain  number  of  electz^^n8. 
Some  of  the  returning  electrons  will  reenter  the  plasma,  cutting 
down  the  net  electron  loss  rate,  but  some  will  be  slightly 
deflected  (a  distance  of  one  gyro  radius)  and  will  miss  the  cusp, 

entering  the  boundary  layer 
(cf.  Pig.  6).  Now,  any  eleo> 
tron  space  charge  in  the 
boundary  layer  decreases 
the  charge  separation  field 
and  allows  the  ion  trajectories 
to  move  out,  thereby  increasing 
Pig.  6  the  ion  loss  rate  (cf.  Pig.  5a; 

the  larger  the  penetration,  the  larger  is  the  distance  d  toward 
the  cusp).  In  addition,  if  this  electrostatic  potential  pene- 
trates inside  the  plasma  through  the  cusp,  it  can  directly  in- 
crease the  loss  rate  of  ions  by  accelerating  them  toward  the  cusp. 

Qualitative  arguments  can  be  given  to  show  that  the  latter 
effect,  i.e.  penetration  of  the  field,  is  small;  this  remains  to  be 
investigated  qxiantitatively.  With  regard  to  the  former  effect,  it  is 
clear  that  the  equilibrium  electron  space  charge  in  the  boundary  layer 
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(I.e.  the  thickness  of  this  layer)  Is  determined  by  a  balance 
between  entering  electrons  which  have  missed  returning  to  the 
cusp  plus  entering  electrons  which  diffuse  from  the  plasma  to 
the  boundary  layer  and  leaving  electrons  which  diffuse  back 
Into  the  plasma  or  out  to  the  walls  of  the  apparatus.   It  Is 
possible  to  conclude  that  a  balance  will  be  reached  In  which 
the  boundary  layer  does  not  thicken  Indefinitely  but  reaches  a 
limiting  thickness  which  Is  not  much  larger  than  the  electron 
gyro  radius.  To  see  this  qualitatively,  we  note  first  that  In 
virtue  of  the  external  electrostatic  field,  electrons  are 
trapped  on  a  magnetic  line  In  the  boundary  layer  whereas  Ions 
are  Immediately  lost  to  the  external  surroiandlngs ;  In  other 
words  the  space  charge  Is  purely  electronic.   Next,  although 
the  electron  density  In  the  plasma  may  be  considerably  larger 
than  that  of  the  electrons  trapped  In  the  boundary  layer,  the 
diffusion  rate  Inward  can  be  considerably  larger  than  the 
diffusion  rate  outward;  (if  not  for  this,  the  boundary  layer 
would  thicken  Indefinitely).  This  happens  because  particle 
deflections  are  predominantly  the  result  of  a  cumulation  of 
small  deflections  rather  than  Isolated  large  deflections.  The 
diffusion  Into  the  plasma  of  a  boundary  layer  electron  Is  such 
a  gradual  affair.  However,  the  outward  diffusion  of  a  plasma 
electron  to  the  boundary  layer  must  be  accomplished  during  a 
single  gyro  period,  while  It  traverses  the  layer;  at  Its  next 
encounter  with  the  layer.  It  has  'forgotten'  any  small  deflec- 
tions It  had  previously  received.   A  rough  quantitative  estimate 
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based  on  these  ideas  yields  an  equilibrixim  boundary  layer  thick- 
ness which  is  not  comsiderably  larger  than  the  electron  gyro 
radius. 

According  to  these  arguments,  the  electric  field  configxiration 
will  be  somewhat  as  shown  in  Fig.  7;  the  lower  figure  is  a  plot 

of  the  electric  potential  as 
a  function  of  radius  along 
the  dashed  line  section  indi- 
cated. 

Making  these  estimates, 
we  obtain  a  formula  for  the 
net  loss  rate, 
(k)  Q  »  10^*^  R  H  ; 

here  Q  is  in  ions  per  second, 
R  is  a  representative  radius 


i ^ *" 

*- boundary     r 

layer 


Pig.   7 

of  the  plasma  in  cm.,  and  H  is  the  magnetic  field  in  gauss. 

This  computation  applies  to  a  geometry  similar  to  Pig.  5c  or  to 

a  single  section  of  Pig.  5d.  The  density  of  the  plasma  is 

implicitly  contained  in  H  since  the  magnetic  field  must  be  such 

as  to  balance  the  plasma  pressure.  The  lifetime  of  the  plasma  is 
given  by 


(5) 


t    -  10'^  H  rVt 


f   is  in  seconds  and  T  is  the  temperature  in  electron  volts, 
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For  example  with  H  «  10,000  guass,  R  -  10  cm. ,  T  -  10  ev.,  we 
have  t  =  0.1  second.   Or,  with  H  =  100,000  gtiass,  R  -  100  cm., 
T  =  10  kev.,  we  have  t  -  0.1  second  again. 

These  nximerical  formulas  are  only  representative;  the 
results  depend  within  a  factor  of  2  or  3  on  which  of  the  various 
cusped  geometries  is  used  (the  above  formulas  are  representative 
of  the  monocusp.  Pig.  4c)  and  by  possibly  even  larger  factors 
depending  on  the  ultimate  effect  of  the  electrostatic  field 
sxirrovinding  the  plasma  on  the  structure  of  the  boundary  layer. 

The  loss  picture  is  entirely  different  in  the  presence  of 
a  crossed  field  (Refs.  (13]  and  [16]).  To  illustrate  this  dif- 
ference, consider  the  following  three  different  particle  loss 
mechanisms: 

(1)  diffusion  across  the  field 

(2)  collision  dominated  mirror  type  losses 
(5)  non-adiabatically  dominated  losses. 

The  first  is  8elf-e3cplanatory  and  is  caused  by  interparticle 
collisions.  The  second  involves  the  concept  of  adiabatic  or 
slow  space  and  time  variation  of  the  magnetic  field  compared  to 
the  gyro  radius  and  frequency.  This  assumption  has  the  conse- 
quence that  a  particle  remembers  its  past  and,  if  once  reflected 
by  a  Biirror,  will  be  in  the  same  situation  at  every  futxire  en- 
counter with  the  mirror,  and  will  be  perpetually  reflected.  This 
ideal  result  is  modified  by  collisions  which  cause  a  particle  to 
forget  its  past.  The  third  loss  mechanism  applies  to  the  same 
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mirror  geometry  when  the  field  changes  appreciably  compared  to 
the  gyro  radius.   In  this  case,  even  without  collisions,  a 
particle  may  be  lost  after  several  traverses  because  of  a  loss 
of  memory.  This  is  exactly  the  sitxiation  In  the  cusped  geometz^^ 
just  analyzed;  the  fact  that  a  particle  is  reflected  at  a  cusp 
yields  no  information  about  the  result  of  a  later  encounter  with 
the  cusp.   We  could  subdivide  case  (?)  further  depending  on 
whether  there  is  a  sharp  interface  (smaller  loss  rate)  or  not 
(larger  loss  rate) — the  difference  is  at  least  one  order  of 
magnitude.  A  non-adlabatically  dominated  geometry  will  usually 
exhibit  losses  which  ai*e  one  or  two  orders  of  magnitude  larger 
than  a  mirror-dominated  geometry,  and  this  in  turn   will  usually 
be  one  or  two  oirders  of  magnitude  higher  than  those  of  a  dif- 
fusion-dominated geometry.  The  total  spread  is  at  least  three 
to  fo\ir  orders  of  magnitude! 

In  general,  all  three  mechanisms  will  be  operative,  but 
frequently  a  single  one  is  dominant.  However,  it  is  even 
possible  for  the  different  mechanisms  to  be  dominant  in  different 
parts  of  the  same  machine.   For  exaB^}le,  consider  a  monocusp 
geometiry  (Pig.  4c)  with  field  and  plasma  intermingled.   A  mirror 
ratio  can  be  assigned  to  ecah  magnetic  line  and  the  nxmerical 
values  will  range  from  say  2  or  5  on  an  outside  magnetic  line 
to  a  very  large  value  on  a  line  which  approaches  the  center  of 
the  apparatus  (where  B  «  0).  The  losses  from  the  outer  layers 


* 
More  precisely,  there  is  only  a  slight  persistence  of  memory, 
Ref.  [13]. 
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will  then  be  collision  dominated  as  in  the  conventional  Mirror 
Machine >  whereas  the  losses  from  the  inner  layers  can  be  con- 
siderably larger  since  they  are  non-adiabatically  dominated 
(see  [13]   for  a  more  complete  discussion). 

The  crucial  feature  of  the  crossed  field  is  that  it  reduces 
the  losses  from  group  (5)  to  either  group  (1)  or  group  (2);  in 
the  case  of  a  sharp  interface  to  group  (1)  and,  more  generally, 
to  gz*oup  (2).  The  explanation  is  roughly  that  in  the  presence 
of  the  crossed  field,  the  resultant  field  is  never  zero,  and 
there  is  a  finite  mirror  z>atio  on  each  magnetic  line. 
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3.   Additional  Considerations 

In  this  section  we  shall  discuss  a  number  of  conpllcating 
features  which  must  be  considered  in  the  analysis  of  an  operating 
cusped  thermonuclear  device. 

First  is  the  possibility  of  an  ion  temperat\ire  which  is 
different  from  the  electron  temperature.   It  is  clear  from  the 
preceding  analysis  that  an  electron  temperature  which  is  lower 
than  the  ion  temperature  reducte  the  net  ion  loss  rate.   As  an 
extreme  example,  if  the  electron  and  ion  velocities  are  equal 
(T  /t_  «  3600),  the  loss  rate  is  smaller  than  that  given  in 
Equation  (4)  by  a  factor  of  about  five.   A  small  effect  in  this 
direction  is  provided  by  Bremstrahlung .   A  more  significant 
mechani&m  is  provided  by  the  selective  heating  of  ions  which 
results  from  some  methods  of  heating  (see  below).  However,  it 
is  unlikely  that  for  relevant  densities  and  temperatures  the 
equipartition  rate  between  ions  and  electrons  will  allow  an 
appreciable  temperature  diffeirence  to  subsist  long  enough  to 
be  useful  as  such.   Nevertheless,  there  is  another  possible 
advantage  in  having  the  electron  temperature  low,  even  for  a 
short  while,  because  with  a  very  small  electron  gyro  radius, 
it  is  more  likely  that  a  plasma  of  the  type  postulated  (with  a 
thin  boundary  layer)  can  be  set  up  experimentally. 

There  are  a  few  aspects  of  the  question  of  heating  which 
are  somewhat  different  from  those  in  other  geometries.   Magnetic 
pumping  as  applied,  for  example,  to  a  Stellerator  geometry 
depends  on  an  irreversible  influx  of  energy  arising  fi»om  an 
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alternating  magnetic  field  variation  near  the  gyro  frequency  of 
the  Ions  or  the  transit  time  throtigh  a  bulge.  A  sliLLlar  situa- 
tion Is  found  In  any  of  the  cusped  geometries,  but.  In  this  case, 
the  frequency  should  be  on  the  order  of  the  transit  frequency  of 
Ions  across  the  apparatus.  The  simplest  model  for  computing 
this  type  of  heating  Is  a  moving  magnetic  wall.  Computations 
based  on  this  model  show  It  to  be  q\ilte  feasible.  However,  this 
computation  Is  too  naive  on  at  least  two  counts.  First  of  all, 
rapid  heating  requires  large  acceleration  of  the  plasma  boxindary 
which  may  lead  to  Instability.  The  restriction  to  stable  values 
of  acceleration  (Equation  (1))  allows  the  moving  wall  to  reach 
on  the  order  of  the  Ion  speed  In  accelerating  over  a  distance 
comparable  to  the  radius  of  the  plasma;  this  allows  the  plasma 
energy   to  Increase  by  a  factor  significantly  greater  than  one 
In  each  cycle,  which  Is  rapid  enough  heating  to  be  practical; 
(any  configuration  which  Is  Interesting  thermonucleairwlse  must 
persist  for  many  Ion  transit  times).  One  particularly  Intriguing 
possibility  for  heating  Is  to  deposit  enoxogh  energy  Into  the 
magnetic  field  to  be  able  to  heat  the  plasma  to  thermonuclear 
energies  In  one  shot  (cf.  the  pinch).  Using  the  same  model,  we 
find  that  the  magnetic  wall  executes  damped  oscillations.  The 
Initial  acceleration  Is,  of  course,  too  high  for  stability. 
However,  within  one  or  two  cycles.  It  drops  to  a  value  on  the 
order  of  the  maximum  stable  value.  This  situation  Is  obviously 
not  one  In  which  the  linear  (I.e.  small  perturbation)  stability 
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analysis  can  be  relied  on  but,  roughly  speaking,  the  possibility 
of  practical  stability.  I.e.  within  a  cycle  or  two,  seems  to  be 
reasonable. 

The  second  major  defect  of  the  moving  wall  model  Is  the 
likelihood  of  creating  shock  waves  In  the  plasma.   In  such  an 
event,  the  moving  wall  merely  accelerates  the  adjacent  plasma, 
and  energy  Is  converted  Into  thermal  energy  of  the  plasma  via 
the  shock.   It  Is  not  definitely  verified  as  yet  that  shocks 
can  exist  \inder  relevant  conditions.   Certainly  a  classical 
shock  whose  width  is  comparable  to  the  mean-free-path  cannot 
exist  since  the  mean- free- path  is  large  compared  to  the  apparatus 
dimensions.  However,  there  are  indications  that  shocks  with 
thickness  on  the  oi?der  of  the  gyro  radius  may  be  observed. 
Asstuning  that  shocks  governed  by  the  gyro  radius  exist,  Hanan 
Rubin  has  shown  that,  for  shocks  of  medlvm  strength,  ions  are 
heated  more  than  electrons  and  under  appropriate  conditions 
even  much  more.   Arguments  have  been  given  by  Stirling  Colgate 
(Ref.  [20])  and  by  Marshall  Rosenbluth  and  Conrad  Longmire 
(Ref.  [19])  to  Indicate  the  opposite  conclusion  for  strong 
shocks.   Insofar  as  all  these  results  are  z>eliable,  they  might 
indicate  a  preference  for  heating  by  a  sequence  of  several 
medium  strength  shocks  rather  than  by  a  single  large  one. 

If  the  situation  is  dominated  by  shock  waves,  the  stability 
analysis  has  to  be  reconsidered.   Instead  of  an  acceleration  of 

♦ 
At  the  present  time  (December  1957)  the  evidence  is  veiry  good 
that  damped  shocks  exist  without  the  intervention  of  particle 
collisions. 
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the  plasma,  we  are  presented  with  an  applied  Impulse  or  sequence 
of  impulses.  This  analysis  has  not  yet  been  done. 

The  single-shot  heating  process  is  basically  the  same  as 
In  the  conventional  pinch,  especially  for  the  geometry  of  Pig.  4a. 
Roughly  speaking,  in  adopting  a  cusped  configviration  for  the 
pinch,  we  achieve  stability  at  the  expense  of  greater  losses. 
The  insertion  of  a  crossed  field  (B„)  increases  stability  in 
both  cases,  and  it  also  cuts  down  losses  in  the  cusped  case. 
We  still  maintain  greater  stability  in  the  cusped  case  than  in 
the  conventional  pinch,  but  we  also,  in  all  probability,  have 
greater  losses  here  as  before.   It  is  difficult  to  Judge  at  the 
present  time  whether  stability  or  particle  losses  is  more  important 
in  obtaining  large  containment  times.   It  is  even  possible  that 
the  relative  importance  could  reverse  when  comparing  scale  models 
and  full  scale  operating  devices. 

The  next  question  concerns  the  extent  to  which  the  velocity 
distribution  is  Maxwellian.  The  cusp  loss  formulas  show  that 
fast  particles  are  lost  more  quickly  than  slow  ones.  This  effect 
will  tend  to  cut  off  the  Maxwell  tail  and  so  reduce  the  number  of 
thermonuclear  reactions.   At  a  given  temperature  in  the  range  of 
10  -  100  kev.,  the  effect  is  much  more  serious  with  DD  than  with 
DT.  Especially  in  preliminary  (low  temperature)  experiments 
when  almost  all  reactions  originate  in  the  tail  of  the  distri- 
bution will  this  effect  be  important.   The  situation  becomes  more 
complicated  \ander  operating  conditions  in  the  presence  of  hot 
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(charged)  reaction  products.   The  primary  effect  on  the  structure 
of  the  boundary  layer  (due  to  the  large  gyro  radius  of  the  hot 
Ions)  Is  to  cut  down  the  loss  rate  of  the  cold  particles.   This 
effect  Is  probably  more  than  enoiigh  to  compensate  for  the 
diminished  Maxwell  tall.   On  the  debit  side,  the  hot  particles 
themselves  are  lost  very  quickly.   First  of  all,  this  makes  vin- 
llkely  the  possibility  of  directly  recovering  the  released 
thermonuclear  energy  by  expansion  of  the  plasma  against  the 
containing  field  since  the  excess  pressure  due  to  the  hot  particles 
Is  small.   Secondly,  this  may  cut  down  on  the  total  thermonuclear 
yield  by  reducing  the  number  of  secondary  reactions  (in  the  case 
of  D-D).   Since  the  computation  of  thermonuclear  yield  Is  fairly 
sensitive  to  the  exact  shape  of  the  Maxwell  tall,  only  rough 
estimates  of  these  factors  can  be  made  at  present.   In  the 
presence  of  a  crossed  field,  the  result  Is  even  more  complex 
because  we  are  faced  with  a  combination  of  non-adlabatlc,  cusp- 
type  losses  which  are  larger  for  hot  particles  and  mirror-type 
losses  which  are  smaller  for  hot  particles.  This  subject  remains 
to  be  investigated. 

The  fact  that  most  of  the  particle  loss  is  in  a  directed 
beam  might  eventually  txim  out  to  be  Important  in  recovering 
this  streaming  energy  directly  as  electrical  energy. 

One  of  the  most  difficult  questions  to  answer  theoretically 
is  the  effect  of  magnetic  field  penetration  into  the  plasma. 
This  factor  is  Important  since  we  do  not  expect  to  easily  construct 
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a  plasma  which  la  practically  field-free.   The  qioantltatlve  loss 
computations  which  have  already  been  done  can  be  definitely 
relied  upon  only  If  the  Internal  gyro  radius  Is  large  compared 
to  the  apparatus  dimensions.   This  is  a  difficult  situation  to 
achieve.   At  the  opposite  extreme,  namely  complete  field  pene- 
tration, we  have  the  picket  fence,  for  which  the  losses  may  be 
much  larger  but  for  which  no  quantitative  results  have  yet  been 
obtained.   The  question  is  how  close  to  the  former  (cusp)  con- 
figuration can  experimental  techniques  bring  us.  With  the 
crossed  field,  there  is  the  same  distinction  between  sharp  and 
gradual  transitions,  but  the  distinction  may  not  be  as  important 
in  this  case,  namely  in  comparing  'good'  with  'better'  rather 
than  'bad' with  'worse'. 

An  Intriguing  feature  is  the  possibility  of  controlling 
the  bovindai*y  layer  thicloiess  by  external  means.  The  particle 
loss  situation  is  improved  if  we  can  draw  off  electrons  from 
the  boundary  layer.   Alternatively,  what  we  wish  to  do  is  lower 
the  potential  of  the  plasma  with  respect  to  the  potential  of 
the  bovindary  layer.  This  goal  would  seem  to  be  possible,  based 
on  the  suggestive  concept  of  high  conductivity  of  a  plasma 
parallel  to  the  field  and  low  conductivity  perpendicular  to  the 
field  (Hall  effect).   Specifically,  Instead  of  the  potential 
distribution  of  Pig.  7,  we  might  be  able  to  obtain  that  of  Pig.  8. 
The  equlpotentlals  would  not  close  around  the  plasma,  but  would 

« 
The  formulas  of  Hartland  Snyder  give  upper  estimates  only. 

-  51  - 


\   /^^^  potential 


Fig.  8 
follow  the  magnetic  lines  and  be  fixed  by  the  extexmal  potentials 
maintained  where  the  magnetic  lines  hit  an  external  conductor. 


-  32  - 


6.   Feaalblllty  Estimates 

An  estimate  of  the  feasibility  of  a  given  geometry  as  a 
thermonuclear  device  can  be  arrived  at  by  balancing  thermonuclear 
power  Input  against  various  losses  to  get  a  positive  net  output. 
As  a  rough  estimate  we  consider  cusp  losses  and  coll  heating 
losses  only,  and  envision  a  steady  state.  The  principal 
omissions  are  radiation  losses  and  the  effect  of  a  duty  cycle. 
We  shall  use  the  particle  loss  estimates  given  in  Section  4  for 
the  raonocusp.  Pig.  4c.  These  are  subject  to  revision,  up  or 
down,  depending  on  the  accuracy  of  the  self -consistent  field 
computation  described  in  Section  4.  More  optimistic  results 
can  be  expected  in  crossed-field  configurations. 

We  can  specify  the  state  of  the  plasma,  once  the  geometry 
is  decided,  by  three  parameters,  e.g.  volume,  pressure  and 
temperature.  Actually  we  choose  R  (a  representative  plasma 
radius  in  cm.),  H  (the  magnetic  field  in  gauss),  and  T  (the 
temperature  in  electron  volts).   For  the  energy  loss  through 
the  cusp,  in  watts,  we  find  (Ref.  [12]) 

(6)  P(,  =  .03  RTH  , 
for  the  ohmlc  losses, 

(7)  Pq  =  10'^  RH^ 

and  for  the  thermonuclear  yield  (charged  particles  only), 
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(8) 


Pj^  -  2  X  lO"-^^  R^'^   ,   T  -  10* 


"   -  2  X  10"^^  r'h*   ,    T  -  10^ 


"   -  5  X  lO"-^"^  rV   ,   T  -  10*^ 


"   -  6  X  lO"^®  rV   ,   T  -  105 


DT 


DD 


Clearly,  P^.  tan  always  by  made  to  dominate  P^  +  P^ 
simply  by  making  R  and  H  large.  Let  us  define  'breakeven' 
by  the  condition  Pj.  «  5(Pp  +  Pr)  (implying  energy  withdrawal 
by  means  of  a  thermal  cycle).   A  few  representative  breakeven 
figures  are  as  follows: 


Reaction 

T 

H 

R 

n(ions/cc) 

r(|i  sec. ) 

ion  transits 

Pc(KW) 

DT 
DD 

10* 
105 

105 
10^ 

22 

40 

3  X  10l« 
3  X  10^8 

5,000 
16,000 

10* 
5  X  10 

6  X  10^ 

io8 

F„     (the  cusp  loss  power)  represents  what  must  be  inserted  to 
generate  the  plasma  and  is  also  comparable  to  the  expected  net 
output  of  the  device.  The  DT  figure  for  power  output  is 
comparable  to  that  of  a  conventional  generating  station,  but  the 
DD  figure  would  seem  to  be  impractical.  The  power  output  scales 
as  1/ yiT  ,  but  H  has  already  been  chosen  as  high  as  is 

feasible  (or  possibly  higher)! 

jf 

This  is  optimistic  in  virtue  of  the  presence  of  a  duty  cycle  but 
pesimistic  in  that  the  energy  supplied  to  create  the  plasma  is 
assumed  to  be  completely  lost. 
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In  the  case  of  a  croased-fleld  configuration,  no  quanti- 
tative particle  loss  results  are  yet  at  hand.   As  a  rough  estimate, 
particle  losses  could  be  computed  by  Interpolation  between 
formulas  representative  of  the  Mirror  Machine  and  the  Stellerator 
geometries.  The  breakeven  figure  for  the  cusped  geometry  would 
then  be  quite  encouraging  since  it  could  be  computed  for  3=1. 
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7.   Possible  Experiments 

In  this  section  we  shall  give  a  brief  description  of  some 
experiments,  not  thermonuclear  In  scope,  which  might  shed  light 
on  the  feasibility  of  constructing  cusped  plasma  configurations. 

One  significant  feature  of  these  geometries  Is  the  sharp 
Interface.  This  can  be  achieved  only  If  the  gyro  radius  is  small 
compared  to  the  size  of  the  apparatus.   As  a  matter  of  fact.  It 
Is  easily  verified  that  the  average  number  of  transits  made  by  a 
particle  before  being  lost  is  directly  proportional  to  the  ratio 
of  plasma  dimension  to  gyro  radius.   As  an  example  of  a  crude 
experiment,  take  the  nominal  Ion  gyro  radius  equal  to  the  plasma 
dimension;  from  the  analysis  of  Section  4,  we  might  expect  an 
optimum  boxindary  layer  thickness  of  about  .05  times  the  plasma 

size.   At  a  plasma  temperature  of  10  volts  and  ion  density  of 

14 
10  ,  we  compute  a  containment  time  of  about  300  microseconds 

which  is  about  4500  electron  transits  or  75  ion  transits.   On 

the  basis  of  free  particle  motion  in  the  vacuum  field  produced 

by  the  external  field  coils  alone  (i.e.  for  a  picket  fence), 

one  would  expect  a  containment  time  of  only  several  ion  transit 

times  or  several  hxindred  electron  transit  times  (this  estimate 

is  tentative  and  must  be  checked  by  numerical  computation  of 

representative  orbits).  This  experiment  should  therefore  allow 

the  two  theories  to  be  differentiated.   With  the  figures  given 

above,  the  field  strength  at  the  plasma  surface  is  about  150 

gauss.  The  radius  of  the  plasma  is  about  5  cm.  The  ion  loss 

-  56  - 


19 
rate  through  the  cusps  is  7.5  x  10   per  second  or  about  10 

amperes.   This  indicates  the  size  of  ion  source  or  plasma  gun 

which  would  be  required  to  create  such  a  plasma. 

It  is  important  to  realize  that,  if  such  a  plasma  is  to 

be  inserted  into  an  already  existing  field  of  appropriate  geometry, 

this  must  be  done  very  quickly.   Otherwise,  the  loss  rate  while 

the  plasma  is  being  created  is  the  larger  value  corresponding 

to  the  picket  fence  mechanism.   A  plasma  gun  rather  than  an  ion 

soxirce  is  therefore  indicated.   A  possible  configuration  is 

indicated  in  Pig.  9- 


source  p 


source 


Pig.  9 
An  alternative  procedure  is  to  first  create  the  plasma  and 
then  apply  the  confining  field.   For  example,  a  shock  can  be 
passed  through  a  tube,  ionizing  the  gas,  after  which  the  magnetic 
field  is  applied.   This  would  seem  to  be  an  efficient  way  of 
constructing  the  desired  sharp  interface.   To  obtain  a  stationary 
plasma,  one  might  utilize  two  shocks,  colliding  head-on  in  the 
center  of  the  tube.   Any  of  the  linear  (i.e.  non-toroidal)  magnetic 
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fields  of  Fig.  4  could  be  applied  to  such  a  plasma  including  the 
ones  with  Internal  (crossed)  fields  and  Internal  conductors 
(holes).   In  addition  to  confining,  the  magnetic  field  could  be 
used  to  compress  the  plasma  and  also  to  remove  it  from  the  walls. 
It  should  not  be  difficult  to  experimentally  determine  the  effect 
of  a  crossed  field  on  containment  and  the  efficiency  of  an 
internal  conductor  for  compression  (see  Section  8). 

Another  method  of  creating  a  plasma  is  by  a  discharge  as 
in  the  pinch.  This  has  the  feature  that  the  same  mechanism 
both  creates  and  then  confines  the  plasma.  The  linear  geometry. 
Fig. 4a,  is  probably  most  appropriate  for  Initial  experiments  of 
this  type.  The  discharge  required  to  produce  this  config\iration 
is  quite  unusual.   It  consists  of  four  parallel  line  discharges, 
alternating  in  direction.  Fig.  10.   It  is  obviously  impossible  to 


four  discharges 


Pig.  10  Pig.  11 

maintain  such  an  array  of  discharges  in  the  absence  of  a  strong 
longitudinal  field,  B^.   With  such  a  field,  it  is  not  clear 
that  the  discharges  would  flow  as  shown  rather  than  the  shorter 
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way,  across  the  magnetic  field,  but  such  a  possibility  does  not 
seem  to  be  ruled  out.  This  experiment  could  probably  be  tried 
rather  easily.   A  preliminary  experiment  might  try  to  run  two 
discharges  separated  by  a  distance  of  several  Ion  gyro  radii. 
The  four  currents  shown  In  Fig.  10  are  eventually  to  be 
surface  currents  on  the  plasma,  (Fig. 11),  and  the  circuits  must 
be  completed  by  Introducing  external  conductors  as  returns.  The 
complete  circuit  might  appear  as  in  Fig.  12;  the  outer  shells 

represent  currents  in  conductors 
and  the  inner  +  signs  represent 
discharges.  This  configuration 
bears  a  strong  resemblance  to 
the  stabilized  (B  )  pinch.   In- 
stead  of  a  single  discharge  and 
concentric  return,  there  are  foxr 
discharges  and  returns.   Now,  the 
Fig.  12  striking  featxare  of  the  stabilized 

pinch  configuration  is  that  ionization  is  produced  by  the  initial 
stages  of  the  same  field  which  later  produces  the  pinch  and  con- 
tainment. Following  this  line.  It  may  be  possible  to  apply  the 
external  (shell)  currents  in  Fig.  12  so  quickly  that  the  con- 
comitant large  electric  field  ionizes  and  initiates  a  discharge 
a  s  shown . 

The  initial  stages  of  this  cusped  B_  pinch  might  take  place 
as  follows.  Upon  application  of  a  potential  (e.g.  by  condenser 
discharge),  the  initial  E  in  the  xinionized  gas  will  be  somewhat 
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(a)  (b) 

Fig.  15 
as  shown  In  Pig.  l^a.  This  Initial  E-fleld  Is  orthogonal  to  the 
B-fleld  that  would  eventually  arise  In  a  vacuum  and  the  magnlcude 
of  E  varies  linearly  along  the  length  of  the  tube  from  zero  at 
one  end  to  the  applied  E  at  the  other;  (It  Is  this  linear 
variation  which  provides  a  curl  E,  consequently  a  c)B/c)t).  Upon 
Ionization,  since  the  current  paths  that  would  be  produced  by 
this  E  Intersect  the  Insulating  walls,  E  Is  Immediately 
reduced  to  zero  Inside  the  plasma  and  Is  concentrated  In  the 
Insulator,  Fig.  13b.  Henceforth  ,  B  grows  only  In  this  Insulator 
and  eventually  builds  up  to  a  value  sufficient  to  push  the  plasma 
away  from  the  wall  as  a  cusped  pinch.   Any  B   that  was  present 
In  the  unionized  gas  Is  trapped  In  the  plasma.  A  current  In  a 
central  conductor  (Pig.  4f )  should  be  Initiated  only  after 
Ionization. 

The  losses  at  the  ends  of  a  cusped  pinch  can  be  cut  down 
considerably  by  using  Insulated  ends  for  the  cylinder.  The  four 
current  paths  on  the  surface  of  the  plasma  will  than  close  In  a 
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thin  layer  of  plasma  near  the  Insulated  end  and  will  pull  away 
to  provide  a  secondai?y  end  pinch  somewhat  as  in  Fig.  Ik.     This 


I 


pinch 


'S^ 

"TT^ 

V 

/ 

/ 

surface   currents 

/ 

t 


Fig.  14 
might  be  useful  in  increasing  the  containment  time  of  a  pre- 
liminary linear  pinch  experiment,  but,  ultimately,  a  toroidal 
geometry  would  probably  have  to  be  resorted  to. 

The  various  modifications  with  internal  conductors  and 
holes  (cf.  Figs.  4f,  4g,  4h)  should  cause  no  difficulty  in 
small  scale  experiments. 
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8.   Comparison  of  Geometries 

In  this  section  we  present  a  semi -quantitative  (frequently 
qualitative)  comparison  between  the  various  cusped  geometries 
and  some  of  the  other,  at  present  further  advanced,  thermonuclear 
devices,  viz.  the  pinch,  the  Stellerator,  and  the  Mirror  Machine. 
In  addition  to  the  major  factors  of  stability  and  loss  rates, 
there  are  questions  Involving  methods  of  heating,  extraction  of 
energy.  Initiation,  ease  of  construction  of  experiments,  and 
many  others  which  can  at  least  be  analyzed  qualitatively  In  the 
absence  of  any  experimental  experience  with  the  cusped  geometries. 

First  consider  geometrical  scaling  of  the  various  cusped 
geometries.   For  the  monocusp  (Fig.  4c),  It  is  clear  that  one 
should  try  to  minimize  the  cusp  perimeter  compared  to  the  volume. 
This  Implies  a  long  thin  shape.   There  Is  a  limit  to  how  far 
one  can  go  in  this  direction  because  the  stability  will  be 
thereby  reduced.   Of  co\arse  this  configxu*atlon  is  always  stable 
against  arbitrarily  small  pertvirbatlons,  but  the  restriction  on 
the  amplitude  of  a  stable  disturbance  becomes  more  severe  as  the 
configxiration  elongates.   Also,  the  amovmt  of  acceleration  that 
can  be  withstood  is  reduced.   A  reasonable  proportion  might  be  a 
length  four  times  the  diameter.   Given  such  an  optimum  shape,  it 
is  clear  that  the  multlcusp  of  Fig.  4d  is  somewhat  better  because 
of  greater  volxirae  per  line  cusp.   On  the  other  hand,  for  a  given 
plasma  volume,  a  long  multlcusp  will  have  a  relatively  small  ratio 
of  volume  to  cusp  length,  so  the  number  of  sections  should  be 
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taken  quite  small,  probably  only  two  or  three,  ending  with  a 
point  cusp  (Pig.  15)- 


Fig.  15 

For  the  linear  geometry.  Fig.  4a,  neglecting  end  losses, 
we  are  led  to  a  short  squat  configuration  if  we  wish  to  minimize 
the  loss  rate  for  a  given  volume  of  plasma.  Exactly  the  opposite 
will  be  tinae  if  we  compute  end  losses  only.  The  optlm\jm  shape 
is  a  length  to  diameter  ratio  about  the  same  as  the  diameter  to 
gyro  radius  ratio.  Or  we  may  cap  the  ends  as  indicated  in 
Fig.  14,  in  which  case  the  optimum  result  will  be  approximately 
cubical.  In  either  case,  the  long,  thin,  open-ended  configuration 
or  the  capped  one,  the  optimization  is  not  quite  as  good  as  in 
the  previous  case.  Fig.  I5. 

Except  in  the  case  of  large  scale  experiments  approaching 
the  thermonuclear  in  scope,  these  differences  in  geometry  would 
probably  be  secondary  to  considerations  of  experimental  simplicity. 

In  devices  which  do  not  include  an  internal  crossed  field, 
there  is  no  point  in  considering  toroidal  configurations  since 
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they  offer  no  advantage  over  Fig.  15«   With  an  Internal  field, 
however,  since  this  field  can  only  be  supplied  in  Pig.  15  by 
inserting  an  internal  conductor  (Fig.  4h),  and  it  is  not  clear 
whether  it  can  be  applied  at  all  with  the  capped  linear  geometry, 
the  most  serious  contender  might  be  the  toroidal  shape  of  Fig.  4b. 
If  insertion  of  an  internal  conductor  should  prove  to  be  feasible, 
the  configuration  of  Fig.  4h  (or  its  equivalent  as  in  Fig.  15) 
would  come  to  the  fore  again.  The  feasibility  of  an  internal 
conductor  may  well  be  answered  differently  for  a  small  scale 
and  for  a  large  scale  experiment,  e.g.  depending  on  the  impor- 
tance of  contaminants  or  of  bombardment  of  the  conductor. 

Finally,  we  remark  that  for  initial,  exploratory  experiments, 
it  may  well  be  possible  to  tolerate  the  open-ended  linear  con- 
figuration. Fig.  4a,  operated  as  a  pinch. 

We  turn  now  to  a  consideration  of  the  scaling  of  tour 
basic  particle  loss  mechanisms.  The  basic  domain  Is  taken  to  be 
a  cylinder  of  length  L  and  radius  R. 

(1)  Unimpeded  end  losses  at  a  rate  ~  rv^  throtigh  an 

2 
area   -vR  .  The  rate  of  loss  of  energy  scales  as 

(2)  Diffusion  losses  perpendicular  to  the  field, 

Pjj  ^  nL/T-'-/^ 

(5)  Mirror  type  losses, 

Pj^  -  n^R^L/T^/^ 
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(4)   Cusp  losses.   For  easy  comparison,  take  the  linear 
geometry  of  Fig.  4a, 

It  Is  illianlnatlng  to  compare  these  energy  loss  rates  with  the 
thermonuclear  Input, 

P^  ^  n^R^L  f(T)   . 

We  have 

^E  *   ^D  *   ^M  *   ^C  *   ^T  "^ 

T^nL  :   l/nR^  j   1  :   T^n^^V  :   T^^^  f  (T)  . 

A  nxamber  of  elementary  conclusions  follow.   Except  In  the  case 
of  mirror  losses,  the  thermonuclear  input  can  always  be  made  to 
dominate  either  by  taking  n  large  or  the  dimensions  large. 
The  factor  n  is  of  most  importance  for  cusp  losses.  The  size 
factor  is  more  important  for  diffusion  and  cusp  losses  than  for 
end  losses.   Also,  devices  in  which  cusp  losses  or  end  losses 
are  significant  would  tend  to  be  operated  at  lower  temperatxu^es 
than  diffusion  or  mirror-loss  dominated  machines. 

In  a  qualitative  analysis  of  this  type  it  is  unnecessairy 
to  consider  radiation  or  ohmic  losses  in  the  field  coils  since 
these  scale  the  same  way  in  all  cases;  quantitatively  there  does 
exist  a  difference  in  that  the  act\ial  numerical  coefficients  may 
be  different. 


.?>..•  ■■^*. 
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The  fundamental  advantages  of  the  cusped  geometries  lie  In 
stability  and  high  p.  The  chief  disadvantages  are  large 
particle  loss  rate  and  Inefficient  use  of  magnetic  field.  To 
be  definite,  let  us  first  compare  the  cusped  linear  pinch  with 
the  conventional  stabilized  pinch.  Problems  of  Initiation  are 
not  radically  different.  With  regard  to  heating,  the  cusped 
pinch  is  superior  in  that  there  is  no  limit  to  the  amo\mt  of 
stable  compression,  and  if  recourse  is  had  to  heating  by  inter- 
mixing of  crossed  fields,  there  is  probably  no  fear  of  losing 
stability.   On  the  other  hand,  the  compressive  efficiency  of 
the  magnetic  field  is  less  because  the  field  is  weaker  at  the 
plasma  than  at  the  coll,  and  this  is  only  partially  compen- 
sated by  the  insertion  of  an  internal  conductor.   Particle 
losses  are  greater  for  the  cusped  pinch,  except  that  for  a 
small  scale  experiment,  the  capped  pinch  would  probably  have 
much  lower  losses  than  an  open-ended  conventional  linear  pinch. 
It  seems  likely  (but  this  is  still  subject  to  verification) 
that  the  cusped  pinch  could  operate  at  somewhat  higher  p  than 
the  stabilized  pinch;  this  is  based  on  an  estimate  that  the 
optimvun  Internal  field  required  to  cut  particle  losses  is  likely 
to  be  lower  than  the  internal  field  required  for  stabilization. 

For  a  preliminary  experiment  in  which  economy  is  unimportant, 
the  cusped  pinch  offers  the  definite  advantage  of  unlimited 
compression  without  endangering  stability  together  with  a 
relatively  low  particle  loss  rate  at  low  temperatures  (small 
compared  to  10  kev)  but  the  disadvantage  of  a  depleted  Maxwell 
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tall  which  might  be  especially  serious  at  low  temperatures. 

With  regard  to  the  Stellerator  and  the  Mirror  Machine, 
the  chief  advantage  of  a  cusped  geometry  would  seem  to  lie  in 
a  considerably  higher  stable  value  of  p.   The  particle  loss 
rate  would  probably  be  higher  for  a  cusped  geometry,  but  this 
is  not  certain  if  crossed  fields  are  resorted  to.   It  is  con- 
ceivable that  the  cusped  configuration  could  operate  in  a  largely 
diffusion-dominated  regime  which  would  bring  the  losses  below 
that  of  the  Mirror  Machine  and  conceivably  even  below  that  of 
the  Stellerator  by  utilizing  a  larger  vol\ame  to  storface  area 
factor  in  a  linear  geometry.  However,  this  is  at  the  present 
time  sheer  conjecture.  The  magnetic  field  utilization  is  likely 
to  be  somewhat  less  efficient  for  a  cusped  configuration  but 
it  might  be  made  more  efficient  by  the  utilization  of  ferro- 
magnetic  cores. 

So  far  as  initiation  by  pre -ionization,  injection,  pinching, 
shock  waves,  or  any  other  method,  there  does  not  seem  to  be  any 
significant  difference  between  cusped  and  other  geometries. 


« 
This  possibility  was  point  out  by  M.  Levlne. 
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